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Abstract 

From perturbation theory, Green's functions are known for providing a simple and 
convenient access to the (complete) spectrum of atoms and ions. Having these functions 
available, they may help carry out perturbation expansions to any order beyond the first 
one. For most realistic potentials, however, the Green's functions need to be calculated 
numerically since an analytic form is known only for free electrons or for their motion in 
a pure Coulomb field. Therefore, in order to facilitate the use of Green's functions also 
for atoms and ions other than the hydrogen-like ions, here we provide an extension to 
the Ratip program which supports the computation of relativistic (one-electron) Green's 
functions in an — arbitrarily given — central-field potential V(r). Different computational 
modes have been implemented to define these effective potentials and to generate the radial 
Green's functions for all bound-state energies i? < 0. In addition, care has been taken 
to provide a user-friendly component of the Ratip package by utilizing features of the 
Fortran 90/95 standard such as data structures, allocatable arrays, or a module-oriented 
design. 



*EmaiI: kovalp@physik.uni-kassel.de, Tel: +49-561-804-4571, Fax: 4-49-561-804-4006. 
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PROGRAM SUMMARY 



Title of program: Xgreens. 
Catalogue number: To be assigned. 

Program obtainable from: CPC Program Library, Queen's University of Belfast, N. Ireland. 
Users may also down-load a tar-file of the program ratipOO.tar from our home page at the 
University of Kassel (http: / /www. physik.uni-kassel.de/fritzsche/programs.htmlj . 

Licensing provisions: None. 

Computer for which the new version has been tested: PC Pentium II, III, IV, Athlon. 
Installations: University of Kassel (Germany). 
Operating systems: SuSE Linux 8.2, SuSE Linux 9.0. 

Program language used in the new version: ANSI standard Fortran 90/95. 

Memory required to execute with typical data: On a standard grid (400 nodes), one central- 
field Green's function requires about 50 kBytes in RAM while approximately 3 MBytes are 
needed if saved as two-dimensional array on some external disc space. 

No. of bits in a word: Real variables of double- and quad-precision are used. 

Peripheral used: Disk for input /output. 

CPU time required to execute test data: 2 minutes on a 450 MHz Pentium III processor. 
Distribution format: compressed tar file. 

Keywords: Central-field Green's function, confiuent hyper geometric function. Coulomb 
Green's function, Kummer function, multi-configuration Dirac-Fock, regular and irregular 
solutions. 

Nature of the physical problem: In atomic perturbation theory. Green's functions may help 
carry out the summation over the complete spectrum of atom and ions, including the (sum- 
mation over the) bound states as well as an integration over the continuum [1]. Analytically, 
however, these functions are known only for free electrons (V(r) = 0) and for electrons in a 
pure Coulomb field (V(r) = —Z/r). For all other choices of the potential, in contrast, the 
Green's functions must be determined numerically. 

Method of solution: Relativistic Green's functions are generated for an arbitrary central-field 
potential V(r) = — Z(r)/r by using a piecewise linear approximation of the effective nuclear 
charge function Z(r) on some grid (i = 1,...,N): Zj(r) = Z^i + Zur. Then, following 
McGuire's algorithm [2], the radial Green's functions are constructed from the (two) linear- 
independent solutions of the homogeneous equation [3]. In the computation of these radial 
functions, the Kummer and Tricomi functions [4] are used extensively. 

Restrictions onto the complexity of the problem: The main restrictions of the program concern 
the shape of the effective nuclear charge Z(r) = — rV(r), i. e. the choice of the potential. 
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and the allowed energies. Apart from obeying the proper boundary conditions for a point -like 
nucleus, namely, Z(r ^ 0) = Z^mc > and Z(r oo) = Znuc — -^electrons ^ 0; the first 
derivative of the charge function Z(r) must be smaller than the (absolute value of the) energy 
of the Green's function, < \E\. 

Unusual features of the program: Xgreens has been designed as a part of the Ratip package 
[5] for the calculation of relativistic atomic transition and ionization properties. In a short 
dialog at the beginning of the execution, the user can specify the choice of the potential 
as well as the energies and the symmetries of the radial Green's functions to be calculated. 
Apart from central-field Green's functions, of course, the Coulomb Green's function [6] can 
also be computed by selecting a constant nuclear charge Z(r) = Z^s- In order to test the 
generated Green's functions, moreover, we compare the two lowest bound-state orbitals which 
are calculated from the Green's functions with those as generated separately for the given 
potential. Like the other components of the Ratip package, Xgreens makes careful use of 
the Fortran 90/95 standard. 
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LONG WRITE-UP 



1 Introduction 

The use of Green's function has a long tradition for solving physical problems, both in classical 
and quantum physics. From perturbation theory, for instance, these functions are known for 
providing a rather simple access to the complete spectrum of a quantum system and, hence, to 
facilitate the computation of perturbation expansions beyond the first order in perturbation 
theory. Applications of the Green's functions can be therefore found not only in atomic and 
molecular physics but also in quantum optics, field theory, solid-state physics, and at various 
place elsewhere. 

In atomic physics, however, the use of Green's function methods was so far mainly restricted 
to describe the motion of electrons in a pure Coulomb field, i. e. to the theory of hydrogen and 
hydrogen-like ions. For these one-electron systems, calculations have been carried out, for 
example, for the two-photon ^ H] and multi-photon ionization the two-photon decay 

0, the second-order contributions to the atomic polarizabilities [£j as well as for determining 
radiative corrections |H1 • 

Less attention, in contrast, has been paid to utilize Green's functions for non-Coulomb 
fields or for describing the properties of many-electron atoms and ions. Unlike to the gen- 
eration of the — bound and free-electron — wave functions to the Schrodinger and Dirac 
equation, for which a number of programs are available now also within the CPC-library 
[HI El ) there is almost no code freely available which helps generate the relativistic central- 
field Green's functions. As known from the literature, however, nonrelativistic central-field 
Green's functions were constructed by McGuire and by Huillier and coworkers and 
were successfully applied for studying the multi-photon ionization of valence-shell electrons 
in alkali atoms. Therefore, in order to facilitate the use of relativistic central-field Green's 
functions for atomic computations, here we describe and provide an extension to the Ratip 
package which calculates these functions (as the solution of the Dirac equation with a 5— like 
inhomogeneity) for an arbitrary central field V(r). 

In the following section, we start with summarizing the basic formulas for the computation 
of relativistic central-field Green's functions. Apart from a brief discussion of the Dirac Hamil- 
tonian, this includes the defining equation for central-field Green's function and the separation 
of the three-dimensional Green's function into radial and angular parts. However, since the 
separation has been discussed in detail elsewhere in the literature we restrict ourselves to 
a short account on that topic and mainly focus on the computation of the radial components 
of the Green's functions. Section |21 later describes the program structure of Xgreens, its 
interactive control and how the code is distributed. Because the Xgreens program is de- 
signed as part of the Ratip package, we have used and modified several modules which were 
published before along with other components of the program. Section |3 explains and displays 
two examples of Xgreens, including (a) a dialog in order to calculate a central-field potential 
from Grasp92 wave functions 9^ and (b) the generation of the radial Green's functions if the 
potential is loaded from an external file. Finally, a short summary and outlook is given in 
section 13 
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2 Theoretical background 



2.1 Relativistic central— field Green's functions 

Most naturally, the relativistic central-field Green's function can be considered as a general- 
ization of the (relativistic) Coulomb Green's function if, in the Dirac Hamiltonian^ 

i?D(r) = -icaV + (3c^ + Vir) , (1) 

the Coulomb potential is replaced by some (arbitrarily given) central-field potential, V^(r) = 
— y — > — As in the nonrelativistic case where the Green's function obeys the defining 
equation ( — £^ ) G^(r, r') = 5{r — r') , the relativistic Green's function is given by a 4 x 4 
matrix which satisfies the inhomogeneous equation 

HB{r) - E - c^) GE{r,r') = 5{r-r')I^, (2) 



with 1 4 being the 4x4 unit-matrix and where, as usual in atomic structure theory, E refers 
to the total energy of the electron but without its rest energy nie c"^ ■ Moreover, since in polar 
coordinates a central-field potential Z(r) in the Hamiltonian does not affect the separation 
of the variables, the central-field Green's function has the same radial-angular representation 
as in the pure Coulomb case ^3] 

^ 1 / gi,^ (r, r') (f) f^L(f') -^gif(r,r')^^™(f)^^L_(f')^ 
GB(r,r') = V— , , (3) 



where riKm(f) = ^^rni'&^'-p) denotes a standard spherical Dirac-spinor and k = ± ( j + 1/2) 
for / = j ± 1/2 is the relativistic angular momentum quantum number; this number carries 
information about both the total angular momentum j as well as the parity (—1)' of the Green's 
function. While the summation over k = ±1, ±2, ... runs over all (non-zero) integers, the 
summation over the magnetic quantum number m = — j, — j + 1, . . . , j is restricted by the 
corresponding total angular momentum. Moreover, the radial part of the central-field Green's 

(g^^ (^ ^') ^')\ 
I in (131) can be treated simply as a 2 x 2 matrix function which 
ginir.r') gjf (r,r')/ 
satisfies the equation 



V 



a Z(r) 



- aE 





K d 




r dr 



d K 
dr r 



a Z(r) 



aE 



^g£« g|f 



a 5{r 



(4) 



with 1 2 now being the 2x2 unit-matrix. Note that in Eq. a refers to Sommerfeld's fine 
structure constant and that, in order to keep the equations similar to those for the Coulomb 
Green's functions ^], we make use of a nuclear charge function Z(r) = — rV(r) to define the 
central-field potential, instead of V(r) explicitly. In the radial-angular representation Q, two 
superscripts T and T' were introduced to denote the individual components in the 2x2 radial 
Green's function matrix. These superscripts take the values T = L or T = S" to refer to 
'^Here and in the following, we use atomic units (me = ft = /A-Keo ~ 1) if not stated otherwise. 



5 



either the large or small component, respectively, when multipHed with a corresponding (2- 
spinor) radial solution of the Dirac Hamiltonian For a pure Coulomb potential, Z(r) = 
Zgfj, an explicit representation of the (four) components g^J' (r, r') of the radial Green's 
function can be found in Refs. [T^ ITU [T3] . 



2.2 Generation of radial Green's functions 



When compared with the Coulomb Green's functions, not much need to be changed for the 
central-field functions in Eqs. Q and ^ except that the nuclear charge Z = Z(r) now depends 
on r and, hence, that analytic solutions to these equations are no longer available. Therefore, 
to find a numerical solution to Eq. let us first mention that this matrix equation just 
describes coupled equations for the two independent pairs (g^^, gj^) ^^'^ Sek) 
radial components. For example, if we consider the first pair (g^^, gj^^) of radial components 
of the Green's function, it has to satisfy the two equations 



a Z(r) 



aE 



IEk 



■ {r, r') + 



d K 
dr r 



.LL 



2 

a 



a Z{r) 



d_ 

dr 

aE 



^SL 

IEk, 



(r, r') = a 6{r — r'), 



lO: (r, r') 



0, 



(5) 
(6) 



and a similar set of equations holds for the second pair (gjf , g^f )• In the following, therefore, 
we will discuss the algorithm for solving Eqs. © and © but need not display the analogue 
formulas for the pair (gjf , g^f )• 

Inserting gj^ (r, r') from Eq. © into ©, we arrive at the second-order inhomogeneous 
differential equation 



d_ 

dr 



d K, 
dr r 



gJ^j^(r,r') 



2 aZ(r) 

- + — + aE 

a r 



+ + Sk^ir,r') = a6{r-r') 



(7) 



for the component g^^(?^, r'). Solution of this equation can be constructed as product of two 
linearly independent solutions jH 



M|^(min(r,r')) • (max(r, r')) 



LL 



(8) 



for the corresponding homogeneous case [cf. Eq. ()1U() below], where (r) denotes a solution 
which is regular at the origin, and (r) a solution regular at infinity. Below, we will obtain 
these functions following the numerical procedure as suggested by McGuire 

For this, let us start with approximating the nuclear charge function Z(r) on some grid 
rj, i = 1, . . . , Zmax in terms of a set of straight lines 

Zj(r) = Zoi + Zur, for ri < r < r^+i , (9) 

and from the homogeneous part of Eq. (O 



d_ 

dr 



d K 
dr r 



2 aZ^ 

— I \- aE 

a r 



+ 



aZj(r) 



+ aE\ ghAr) = 



(10) 
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as given for the i-th interval [r,, rj+i] of the grid. In this Eq., we can drop the second 
argument r' in the component g^^, since it now appears only as a parameter, and replace it 
by the superscript i to denote the particular piece of the grid for which we want the solution. 
From the approximation @, moreover, we see that the nuclear charge function Zj(r) within 
the i—th interval gives rise to a pure Coulomb potential Zoi/r and a constant which is simply 
added to the energy: E — > Z^i + E. We therefore find that, within the given interval, the 
regular and irregular solutions at the origin in ((HJ can both be written as linear combinations 
of the corresponding solutions for a Coulomb field 

Mh^ir) = AiM^^C— (r) + /,2W^,^°"(r) ^^^^ 

W^Jr) = 5oM^,^°"(r) + g,,2W^f"(r) (12) 

with constants /j,2, 9i,i, 91,2} which still need to be determined. For the Coulomb po- 

tential, the functions M^^ °" °™ (r) and W^^ °" °™ (r) are known analytically and given by 

M^^^°"(r) = r^' e-«>^ [U M(-i, + 1, Is^ + 1, 2q,r) 

+ {k- Zoi/q,) M{-U, 2si + 1, 2q,r)] , (13) 

W^^^°"(r) = e-"'^ [{k + Zo./<7i) U(-t, + 1, 2si + 1, 2g,r) 

+ U(-t„2s, + l,2(7ir)] , (14) 

where M(a, 6, r) and U(a, 6, r) denote the Kummer and Tricomi functions |161 117j . respec- 
tively, and where the quantities Sj, U and qi are given by 



aZoi{{E + Zu)a^ + l) 



VI - {{E + Ziija^ + 1) 



Qi = V-iE + Zu){iE + Zu)a^ + 2) . (16) 

The set of constants /i,2, 5j,2; ^ = 1, •••) ^max} can be determined from the fact 

that the two functions M^^(r) and Wj;^(r) in ansatz Q as well as their derivatives need 
to be continuous in r, and that they behave regularly at the origin or at infinity, respectively. 
The constraint of being regular at the origin, for instance, requires the coefficients /j=i,i = 1 
and /i=i,2 = and can be used together with the continuity M^^(r) and M^^-^(r) , 

h,iM^^{r,) + n,,W^'^{n) = M^+i'C(r,) + /m,2T^j^:'''^(r.), (17) 

hiM^'fir.) + n,2W;;f(r,) = /.+i,iM;;^^'^(r.) + /.+i,2t^il+^'^(r.). (18) 

[where the superscript C here refers to the Coulomb functions in Eqs. (|13|) and (|14)) ] in order 
to determine all the coefficients fij up to a normalization constant. A similar recurrence 
procedure also applies to the coefficients gi,j, but by starting from 'infinity', that is with 
5imax,i = 5'imax,2 = 1) ^ud by golug bttckwards in the index i towards the origin. 

To determine finally the normalization of the radial component g^^(r, r') , e. g. of the 
coefficients {fij, gi,j}, we may return to Eq. Q and re-write it in the form 

ho - e) gi^(r, /) = 6{r-r'). (19) 
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Taking the integral over latter equation, we see that 



i-r'+e ^ rr'+e 

/ {ho - E) g^Ki^^r') dr = / 6{r-r')dr 

Jr'—e Jr'—e 



should be valid for e ^ +0 and for any r' and, hence, that — by using Eq. ((71) and carrying 
out some algebraic manipulations — the derivative 



§;sti^^ r') 



r = r'+e / r. rj i i\ 

' 1 a L\r ) 



a 



a - + ^ + aE (20) 



'jumps' at r = r' . We can use the right hand side of Eq. (|20j) together with the derivative 
of Eq. © to determine the normalization constants Cf and Cg with which the coefficients /jj- 
and gj.j need to be multiplied in order to obtain the normalized radial component g^^(?', t') . 

Having constructed g^^ (r, r') piecewise as solution to Eq. ((71), we may obtain the second 
component gj^ (r, r') of this pair simply from Eq. 



->Ek 



= > (21) 



2 aZ(r) 
a r 



where the derivatives of the Kummer and Tricomi functions in expressions (|13() and (|14|) can 
be calculated by means of standard formulae ^7] 

M'(a, h,z) = '^ M(a + 1, 6 + 1, z) , U'(a, h, z) = -aU(a + 1, 6 + 1, z) . 

In practise, both functions M(a, 6, z) and U(o, 6, z) are required only for real arguments o, 6, 
and z but need to be calculated with a rather sophisticated algorithm in order to ensure 
numerical stability and to provide sufficiently accurate results. A similar numerical procedure 
has to be carry out also for the second pair (gjf , g^f ) of radial components. This gives 
rise of course to another set of coefficients Qij} and, finally, to the full radial Green's 

function from Eq. Q . The 4 i max coefficients for the (piecewise) regular and irregular solutions 
in ansatz (|llH12j) and (jHI certainly provide — together with a few numerical procedures — 
the most compact representation of the radial central-field Green's functions. For the further 
computation of matrix elements and atomic properties, however, these radial components are 
usually represented (and stored) at some (2-dimensional) grid in r and r' as we will discuss 
below. 



2.3 Tests on the accuracy of the radial Green's functions 

To make further use of the Green's functions in applications, it is necessary to have a simple 
test on their numerical accuracy. For the radial-spherical representation of these functions 
as defined in Eq. ©, such a test is easily constructed since the Green's function contains the 
information about the complete spectrum of the Hamiltonian and, hence, about any of its 
eigenstates V'n (r) . Making use of the well-known expansion of the Green's function 
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in terms of the eigenstates of the Hamiltonian, the relation 

Mr') = {En-E) j dvMv)GE{v, r') 



(22) 



can be easily derived from the orthogonality of the eigenfunctions and holds for all energies 
E 7^ En- With the representation Q in mind, the relation (|22j) can be written also in terms 
of the radial wave and Green's functions 



Ann (r') 



{Enn — E) 



dr (Pn.(r), Q„,(r)) 



(23) 



where P„k (r) and Q„k ('') are the large and small components of the (Dirac) radial 2-spinors. 
As indicated by the tilde on the left-hand-side of Eq. (|23|) . therefore, this relation may serve 
as a test on the accuracy of the Green's functions if, for instance, the radial components from 
both side of the equation are compared with each other or if some proper overlap integral is 
calculated. 

To test on the precision of the Green's functions in Xgreens, we generate the radial 
bound-state wave function PnK(r) and Qlnni'i') as solution of the radial Dirac equation 



a Z(r) 



aE 



d K 
dr r 



r + 



2 aZ(r) 



d_ 

dr 

aE 



QnK(^) 











(24) 
(25) 



with the same potential as applied before and by making use of the program by Salvat 

etal UOj. This solver has been embedded in our code and is utilized in order to calculate the 
relativistic components. With these radial functions, we then compute the two integrals in 
H23|) for obtaining Pnn (f) and Qnn (f), respectively, as a function of r. As default in Xgreens, 
we use the overlap integral 

nK(^) PfiK 

(r)+Q„.(r)Q„,(r)) dr (26) 

for the two lowest principal quantum numbers n of the (given) symmetry n, together with the 
corresponding normalization integral 



( P„K(r)P„.(r) + Q„,(r)Q,,(r) ) dr , 



(27) 



as a numerical measure on the accuracy of the generated Green's functions. These integrals 
are displayed explicitly by the program (on demand). For a standard (logarithmic) grid with 
about 300 nodes, the overlap integral (|26() is typically within the range 10~^ . . . 10~^. For 
such a grid, the accuracy is limited in the computations by the linear interpolation of the wave 
and Green's functions. The accuracy can be increased however for a larger number imax of 
grid points as shown in section [l] 



3 Program structure 
3.1 The RATIP package 

Similar as Grasp92 [2] was designed for generating the wave functions within the multi- 
configuration Dirac-Fock (MCDF) model, the Ratip package |19| I2()| I21j is organized as a 
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suite of program components to calculate a variety of (relativistic) atomic transition and 
ionization properties. Among other features, the components of Ratip support for instance 
investigations on the autoionization of atoms and ions, the interaction with the radiation 
field, the parametrization of angular distributions in the emission of electrons and photons, 
or the analysis of interference effects between radiative and non-radiative processes. In quite 
different case studies, Ratip has helped analyze and interprete a large number of spectra and 
experiments. In order to provide efficient tools for atomic computations, it also incorporates 
a number of further components (beside of the main program components for calculating 
certain physical properties) which facilitate the transformation of wave functions between 
different coupling schemes or the generation of continuum orbitals and angular coefficients. 
With the development of the Xgreens program, we now provide an additional component to 
generate the relativistic central-field Green's functions within the Ratip environment. This 
is rather independent of their later 'use' where different properties such as the two-photon 
ionization and decay or various polarizabilities of atoms and ions might be calculated with 
the help of these functions. Owing to the large number of possible applications, however, a 
support of certain properties will largely depends on the requests by the users and on our 
further experience with the code. 

Not much more need to be said about Ratip 's overall structure. A detailed account on its 
present capabilities has been given previously p2j. With regard to the further development 
of Ratip, we just note that our main concern now pertains to a long life-cycle of the code 
and to an object-oriented design within the framework of Fortran 90/95. With the present 
set-up of the Xgreens component, we continue our effort for providing an atomic code which 
is prepared for future applications in dealing with open-shell atoms and ions. 

3.2 Data structures and program execution 

Certainly, the major purpose of Xgreens is the computation of (one-electron) radial Green's 
functions for some — specified or externally given — central-field potential, which can later 
be utilized for calculating radial matrix elements of the type 

U^'^r^ (fe. A; E- ~k,A) = I gl'{r) ]^{kT) gH (r, r') y^{kr') g^'^(r') dr dr' , (28) 

where jA(A;r) denotes a spherical Bessel function, g^J (r, r') one of the components of the 
radial Green's function in Q, and where a and (5 refer to the radial wave functions of some 
bound or free-electron state, respectively. Matrix elements of this type play a key role, for 
instance, for studying two-photon ionization (bound-free) or decay (bound-bound) processes. 
Since the (four) components of the radial Green's functions depend on two radial coordinates, 
r and r', a rather large amount of data need usually to be calculated by means of Xgreens. 
These data are finally stored in an external (.rgf) radial Green's function file as will be 
discussed in the next subsection. 

To generate and utilize the Green's functions, two derived data types TGreens_single_rgf 
and TGreensjrgf have been introduced as shown in Figure ^ Using these data types, a single 
Green's function is kept internally in terms of its expansion coefficients from Eqs. and 
H12|) for the large-large and a similar set of coefficients for the small-small component on the 
given grid. Apart from the particular representation of the grid, these data structures contain 
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type : : TGreens_single_rgf 

! Data structure type to keep a single Green's function 

real (dp) :: energy 

integer : : kappa 

! Coefficients for the regular and irregular solutions 
real(dp), dimension( : , : ) , allocatable :: fLLl, fLL2, fSSl, fSS2 
! Values of the regular and irregular solutions and their derivatives 
real(dp), dimension( : ) , allocatable :: mL, wL, mS, wS, & 

mLp, wLp, mSp, wSp 

end type TGreens_single_rgf 

! 

! 

type : : TGreens_rgf 

! Data structure type to save the Green's fmictions 

type (TGreens_single_rgf ) , dimension( : ) , allocatable :: gf 
integer :: interpolation_mode , mtp 

real(dp), dimensionC : ) , allocatable :: yp, zO, zl 
end type TGreens_rgf 

Figure 1: The derived data structures TGreens_single_rgf and TGreens_rgf to generate, store, and 
interpolate a central-field Green's function. 

the energy and symmetry of the radial Green's function, the nuclear charge function(s) as well 
as information about the mode of interpolation of the radial components between the grid 
points and the maximum tabulation point mtp. In Xgreens, a variable of type(TGreens_rgf) 
is used in order to store the information about all the requested functions, and a few additional 
procedures are provided to calculate from these structures the values of the radial components 
for any set of arguments. 

The basic steps in the execution of the program are very similar as for other components 
of the Ratip package. At the beginning, all the necessary input data are read in by an 
(interactive) dialog. The requested Green's functions are then calculated in turn and tested 
for their accuracy. The (main) output of the program is a formatted Ascil file which provides 
an interface for further applications, either within the Ratip environment or as worked out 
by the user. 

3.3 Interactive control and output of the program 

Like the other components of the Ratip package, Xgreens is controlled by a dialog at the 
beginning of the execution. In this dialog, all information need to be specified for obtaining the 
central-field potential as well as about the number and symmetry of the radial Green's func- 
tions to be generated. Owing to various choices for defining the spherical potential, however, 
slightly different dialogs may occur during the execution. An example is shown in Figure |21 
where first a Hartree-plus-statistical-exchange (HX) potential due to Cowan |23| I24j is gen- 
erated from the ground-state wave functions of atomic gold, and then two radial Green's 
functions are calculated within this potential. 

Apart from the nuclear charge, given in terms of a Grasp92 .iso isotope data file, the 
dialog first prompts for specifying (or confirming) some basic parameters such as the speed of 
light, the grid parameters as well as for the energy units, in which the further input and output 
is done by the program. This is followed by a prompt for determining the central-field potential 
for the radial Green's functions. At present, Xgreens supports four models including a pure 
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Coulomb potential (for the calculation of Coulomb Green's functions) as well two potentials 
(HFS and HX) to incorporate also the exchange interaction in some approximate form. In the 
Hartree model, in contrast, only the direct potential of some atomic level is applied in the 
later computation of the Green's functions. In all these cases, however, the wave functions 
of the selected state have to be specified in terms of the corresponding Grasp92 files for 
the configuration state functions (CSF), mixing coefficients as well as the radial orbitals. If 
computed internally, the potential can be written also to disc and used in further applications. 
Beside of the internal set-up of the potential, in addition, a user-defined potential can be 
given also explicitly and read in by the program. 

Although the nuclear charge is specified by means of a . iso file (in line with all other 
components of Ratip), the Green's functions are always generated for a point-like nucleus in 
order to 'utilize' their regular behaviour of the radial components at the origin [cf. section 1221 • 
For each radial Green's function, then the energy and the one-particle symmetry k need to 
be specified. The functions are finally written to the .rgf radial Green's function file whose 
name has also to be given at the end of the dialog. In this file, the radial Green's functions 
are given at the grid (in r and r') as specified originally. Moreover, to test the accuracy of the 
Green's functions, a number of overlap and normalization integrals are calculated on request 
in order to 'compare' a few low-lying orbitals, as obtained from the Green's functions by an 
integration over r', with independently generated functions from Salvat's program ^Hj- In 
average, the computation of a single Green's function requires about 2 min on a 450 MHz 
Pentium PC and approximately the same time for the 'test' integrals. In the latter case, the 
rather large demand on CPU time arises mainly from the 2imax radial integrals in r'— space, 
which need first to be performed in order to obtain the radial orbital functions from the given 
Green's function components. 

3.4 Distribution and installation 

A program of Ratip 's size cannot be maintained over a longer period without that certain 
changes and the adaptation of the code to recent developments (in either the hardware or 
the operating systems) become necessary from time to time. For the distribution of the code, 
therefore, we follow our previous style in that the Ratip package is provided as whole. The 
main emphasis with the present extension of the program, however, is placed on the design 
and implementation of the two (new) modules rabs_greens . f 90 and rabs_utilities_2 . f 90 
which contain the code for the generation and the test of the radial Green's functions. In 
addition, the four modules (rabs_special_f unctions . f 90, rabs_error_control . f 90, 
rabs_greens_external . f 90, and rabs_cpc_salvat_1995 . f ) became necessary and had to be 
appended to the code in order to support an accurate computation of the hypergeometric and 
a few related functions from mathematical physics as well as for providing Salvat's solver ^U] 
for the generation of the radial wave functions (P(?'), Q(?')) within a given potential. 

In the present version, the Ratip package now contains the source code for the 8 components 
Anco, Cesd, Lsj, Greens, Rcfp, Relci, Reos as weh as the Utilities for performing a 
number of small but frequently occurring tasks. Together with the six additional modules from 
above, the overall program therefore comprises about 40.000 lines of code, separated into 22 
modules (apart from the main program components and three libraries). All these source files 
are provided in the Ratip root directory. As before, there is one makefile for each individual 
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component from which the corresponding executable can be obtained simply by typing the 
command make -f make-component, that is make -f make-greens in the present case. For 
most components of the Ratip package, in addition, we provide a test suite in a subdirectory 
test-component of the root where component refers to the names above. For example, the 
directory test-greens contains all the necessary files to run the sample calculations which we 
will discuss below. 

Before, however, the makefiles can be 'utilized', a number of global variables need to be 
specified for the compilation and linkage of the program. In a Linux or Unix environment, 
this is achieved by modifying (and sourceing) the script file make-environment which saves 
the user from adopting each makefile independently. In fact, the script make-environment 
just contains a very few lines for specifying the local compiler, the options for the compiler, 
as well as the local paths for the libraries. Under Windows, in contrast, not much help need 
to be given to the user since most compiler nowadays provide the feature to define 'projects'. 
In this case, it is recommended to read off the required modules and source files from the 
makefile of the program component and to 'declare' them directly to the project as associated 
with the given component. 

In the past year, the Xgreens program has been applied under several LiNUX systems 
using the Lahey Fortran 95 compiler. The Ratip program as a whole has been found also 
portable rather easily to other platforms such as IBM RS/6000, Sun OS, or to the PC world. 
The file Read. me file in the ratip root directory contains further details for the installation. 
Overall, however, it is expected that it will not be difficult to compile the program also under 
other operating systems. 

4 Test calculations 

To illustrate the use of the Xgreens program, we briefly discuss and display two examples. 
They refer to the generation of the radial Green's function for a few selected (one-electron) 
symmetries in atomic gold [Z = 79) and demonstrate how the accuracy of these functions 
can be increased by a proper choice of the radial grid. However, before the radial components 
can be generated, we need first to specify the potential as shown in Figure [21 For the present 
examples, we have chosen a HX (Hartree-plus-statistical-exchange) due to the work of Cowan 
[231 124j . starting from the ground-state wave functions of neutral gold. As usual for Grasp92 
and Ratip, these wave functions have to be provided in terms of the .csl configuration state 
list and the .mix configuration mixing files as well as the .out radial orbital file. All of these 
input files are provided also in the test-greens/ subdirectory of the Ratip root directory 
and, thus, can be used for the present tests. In Figure 121 moreover, the potential is saved to 
the central-field potential file z79-au-hx.pot and later re-utilized in Figure IHl The output of 
the two examples are shown in the Test Run Output below. 

Like for the wave functions, the 'quality' of the generated Green's functions becomes only 
fully apparent, if they are used for calculating observables which can be compared to experi- 
ment. Although calculation of observables is not the aim of the present work, it is still possible 
to test the accuracy of the generated Green's functions by re-calculating the radial orbitals (as 
discussed in section IT^ and by comparing these orbitals with those as obtained from the inte- 
gration of the Dirac equation for the given central field. This 'comparison' is done in Xgreens 
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XGREENS: Calculation of relativistic central-field Green's functions 
with energies E < (Fortran 95 version) ; 
(C) Copyright by P Koval and S Fritzsche, Kassel (2004). 

Enter a file name for the greens. sum file: 
z79-au-pot . sum 

Enter the name of the isotope data file: 
z79.iso 

loading isotope data file . . . 

... load complete; 
The physical speed of light in atomic units is 137.0359895000000 ; 

revise this value ? 

n 

The parameters of the radial grid 1 (first grid) are 



rnt = 2 . 177968408335618E-04 

h = 6.250000000000000E-02 

hp = O.OOOOOOOOOOOOOOOE+00 

n = 390 ; 



revise these values ? 

n 

Which units are to be used to enter and to print the energies: 
A : Angstrom; 

eV : electron volts; 

Hartree : Hartree atomic units; 
Hz : Hertz; 

Kayser : [cm**(-l)]; 

eV 

Determine type of central-field potential to generate: 
Load : Load a potential from a .pot file 
Coulomb : (pure) Coulomb potential 
Hartree : Direct potential for a given ASF 
HFS : Hartree-Fock-Slater method 

HX : Hartree-plus-statistical-eXchange (Cowan 1967) 

HX 

Enter the name of the .csl (configuration states list) file: 
au. csl 

Loading configuration symmetry list file . . . 

There are 22 relativistic subshells; 

there are 3 relativistic CSFs; 

. . . load complete . 
Enter the name of the GRASP92 mixing coefficient file: 
z79-au.mix 
Loading mixing coefficients file . . . 

... load complete; 
Enter the name of the rwf data file: 
z79-au.rwf 

. . . load complete ; 

Summary of all ASF to facilitate the selection of proper bound-state levels 
i Level J"P Energy (a.u.) Energy (eV) 



1 1 1/2 + -1.90397E+04 -5.18096E+05 
Enter ASF level number for which the potential should be calculated: 

1 

Enter the subshell (e.g. Is, 2p-, ...) which is specific to Cowan's HX method: 

Is 

Write the generated potential to disc ? 

y 

Enter the name of .pot central-field potential file: 
z79-au-hx . pot 

Radial Green's functions will be calculated within the range <= r, r' <= 4.797 a.u. 
with the boundaries of the nuclear charge function Z(0) = 79.00 and Z(r_max) = 1.951 
In the given potential the one-electron energies (in eV) are: 
Is: -8.131E+04 

2s: -1.453E+04 2p-: -1.402E+04 2p: -1.214E+04 

3s: -3.479E+03 3p-: -3.232E+03 3p: -2.814E+03 3d-: -2.382E+03 3d: -2.289E+03 

Enter in turn the energies E (in eV) and symmetries of the radial Green's functions to be generated; 
Enter another energy and symmetry (<cr> if done) : 
-10000 s 

Enter another energy and symmetry (<cr> if done) : 
-15000 d- 

Enter another energy and symmetry (<cr> if done) : 

Check the accuracy of the Green's functions after generation ? 

y 

Enter the name of .rgf file to save the Green's function: 
z79-au-pot . rgf 



Figure 2: Interactive dialog for calculating the effective potential and for generating the radial Green's 
functions. 
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The parameters of the radial grid 1 (first grid) are 
rnt = 2.177968408335618E-04 ; 
h = 6.250000000000000E-02 ; 
hp = O.OOOOOOOOOOOOOOOE+00 ; 
n = 390 ; 
revise these values ? 

y 

Enter rnt: 
2 . 177968408335618E-04 

Enter h: 
3.125E-2 

enter hp: 



enter n: 
780 

Which units are to be used to enter and to print the energies: 



Angstrom; 
electron volts; 
Hartree atomic units; 
Hertz; 
[cm**(-l)] ; 



A 

eV 

Hartree 
Hz 

Kayser 
eV 

Determine type of central-field potential to generate: 
Load : Load a potential from a .pot file 
Coulomb : (pure) Coulomb potential 
Hartree : Direct potential for a given ASF 
HPS : Hartree-Fock-Slater method 

HX : Hartree-plus-statistical-eXchange (Cowan 1967) 

Load 

Enter the name of a .pot file with a central-field potential 
z79-au-hx.f Id 

Radial Green's functions will be calculated within the range <= r, r' <= 5.269 a.u. 
with the boundaries of the nuclear charge function Z(0) = 79.00 and Z(r_max) = 1.000 
In the given potential the one-electron energies (in eV) are: 
Is: -8.131E+04 

2s: -1.453E+04 2p-: -1.402E+04 2p: -1.214E+04 

3s: -3.479E+03 3p-: -3.232E+03 3p: -2.814E+03 3d-: -2.382E+03 3d: -2.290E+03 
Enter in turn the energies E (in eV) and symmetries of the radial Green's fimctions to be generated; 



Enter another 


energy 


and 


symmetry 


«cr> 


if 


done) 


-1000 s 














Enter another 


energy 


and 


symmetry 


«cr> 


if 


done) 


-10000 s 














Enter another 


energy 


and 


symmetry 


«cr> 


if 


done) 


-15000 p 














Enter another 


energy 


and 


symmetry 


«cr> 


if 


done) 


-10000 p- 














Enter another 


energy 


and 


symmetry 


(<cr> 


if 


done) 


-15000 d- 














Enter another 


energy 


and 


symmetry 


(<cr> 


if 


done) 


-10000 d 














Enter another 


energy 


and 


symmetry 


(<cr> 


if 


done) 



Check the accuracy of the Green's functions after generation ? 

y 

Enter the name of .rgf file to save the Green's function: 
z79-au-rgf . rgf 



Figure 3: Dialog for calculating the radial Green's functions with an increased number of radial points. 
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on request by evaluating for the two lowest principal quantum numbers n the (radial) overlap 
integrals of the re-calculated functions with the solutions by Salvat's program. In addition, 
the normalization integrals are also displayed as seen from the Test Run Output. Except 
of a few cases, the deviation of these integrals from 1.0 are typically well below 0.1 % on the 
standard grid with about 300 mesh points in r and r'. For these deviations, there are two 
source of numerical 'inaccuracy' which arise from the test procedure and are not related to the 
generation of the radial Green's functions. They are caused by the fact that by using Salvat's 
solver ) a third-order spline is used to represent the potential at intermediate points in r in 
contrast to the linear representation Q for the generation of the Green's functions. Moreover, 
the integration of the radial integrals uses the standard Grasp92 procedure and, hence, is not 
adapted as well to Salvat's solutions. 

The accuracy of the Green's functions (and the test integrals) can be improved however 
if the number of grid is enlarged, i. e. the spacing between the mesh points reduced. In our 
second example, cf. Figure |31 we therefore adopt a grid with approximately twice the number 
of mesh points. Note that this requires four times as much storage as the radial components 
are functions in r and r'. In this test case, moreover, the radial Green's functions are generated 
for four different symmetries s, Pi/2) P3/21 ^3/21 When the Test Run Output from 

this example is compared with those from example 1, we see that the accuracy in the overlap 
integrals in increased by about a factor of 5. For this example, the generation and test took 
about 30 minutes in total on a 450 MHz Pentium III. 

After the termination of Xgreens, the radial Green's functions file (e. g. z79-au-rgf .rgf 
from example 2) contains all the Green's function components which were generated during 
the execution. Since this is a formatted file, it can easily be read and manipulated by any 
text editor. This (radial Green's function) file stores the individual components g^J' {r, r') of 
the Green's functions as 2-dimensional arrays within a simple file structure. An (internal) file 
signature # DCFGF in the first line is followed by the mode of interpolation and the number of 
Green's functions which are provided by the file. Then, each Green's function is specified (in 
turn), starting with the energy and symmetry of the function and followed by a table of six 
rows where the first two rows refer to the coordinates r and r' (in atomic units) and and the 
other rows to the four component functions ^E^{T,r')^ SekC'^' ''O) §J'k('^;'"') g§^{r,r') in 
this particular order. Although this format is convenient for the later use of these functions, 
it may lead to rather large data files, especially if some larger number of Greens' functions 
need to be generated. In practise, however, not much problems are expected with this file 
structure because disc storage became cheap recently and because, for most applications, we 
expect the Green's functions to be generated on demand by making use of the corresponding 
modules instead of keeping them in external files. Typically, the radial wave and Green's 
functions occur as part of matrix elements and, thus, first require an additional integration 
(over r and/or r') before any observable quantity is obtained. 

The format of the radial Green's function (.rgf) file can be used easily in order to plot 
the various radial components. In Figure |2 therefore, we display the nuclear charge function 
Z(r) and the two radial Green's function components g^^{ro,r) and g§^{rQ,r) as function 
of (second) argument r and taken for tq = 0.204076 a. u. The radial Green's functions are 
calculated in the ground-state potential of atomic gold, as applied in example 2, for the energy 
E = 1000 eV. As seen from Figure HJ the large-large component g^^ is continuous at r = r' 
while the small-large component gj^ jumps owing to Eqs. (|2()|) and (|21j) at this position 
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Figure 4: Nuclear charge Z = Z(7-) and the radial Green's function components g^^(ro,r) and 
gJ^(ro,r) as function of r and for tq — 0.204076 a.u. The radial Green's functions with s— symmetry 
(k = —1) are calculated in the ground-state potential of atomic gold for the energy E — 1000 eV. 



in space. Both components, moreover, behave 'regularly' at the origin as was constructed 
exphcity by Eq. 

5 Outlook. Applications of central— field Green's functions 

Apart from fundamental interest in having the central-field Green's functions available for 
'relativistic' electrons, these functions are also useful for the perturbative treatment of atoms 
and ions. As discussed in subsection [231 namely, the Green's function provide a simple access 
to the complete spectrum of a quantum system and, hence, can be utilized for carrying out 
the summation over all the (one-particle) states of the spectrum as required in second- and 
higher-order perturbation theory. However, since we only provide the one-electron Green's 
functions, they can be applied for just those processes where the 'perturbations' are described 
in terms of one-particle operators. Perhaps, the most-studied perturbation of this type is the 
interaction of atoms and ions with the radiation field, which — within a sufficiently strong 
field — may lead not only to atomic photoionization but (in second order) also to the two- 
photon ionization and decay of atoms and ions and to several other processes. During the last 
two years, we utilized the central-field Green's functions mainly for exploring the two-photon 
ionization for the helium-like ions and for the inner-shell electrons from atomic neon and 
argon. 

Apart from the one-particle Green's functions (as discussed in this work), there are a 
number of further processes such as the double Auger decay, for which the two~particle Green's 
functions G£;(ri,r2; r'i,r2) are required in order to calculate the cross sections, decay rates 
and/or angular distributions. These two-particle functions are solutions to a defining equation 
similar to Eq. but where the (one-particle) Dirac Hamiltonian Hj) is replaced by the two- 
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particle operator -ffo,! + -f^D,2 + Vu, including the electron-electron interaction V12 explicitly. 
Even by making use of a proper radial-angular decomposition of such Green's functions, their 
radial part g^^ ^2(^1' '^2; ''ii ^2) would depend then on four radial variables and an overall rank 
k, similar as known from the tensorial decomposition of the electron-electron interaction j25j . 
Therefore, an internal representation and generation of these (radial) functions appear rather 
infeasible in practise. An alternative to these two-particle Green's functions is given by the 
(so-called) modified pair functions ^£;,a(ri,r2) which satisfy the equation 

[Hb,i + Hd,2 + V12 - e} $i^,«(ri,r2) = W \^^{PJM)) , 

where is a one- or two-particle operator which depends on the physical task to be solved 
and \^a{PJM)) is one of the (two-electron) solutions of the corresponding homogeneous 
equation. These modified pair functions have the advantage that they only depend on two 
radial variables similar as the one-particle Green's functions. On the other hand, of course, 
these functions now depend on some (initial) two-particle state \^a{PJM)) as well as on the 
particular choice of the interaction operator W and, hence, are less general when compared 
with the Green's function from above. We currently investigate the possibilities for generating 
also such modified pair functions within the framework of Ratip and for their (later) use in 
the description of the double Auger decay. But already the one-particle Green's functions 
from this work enables one to explore a number of atomic properties which have not been 
studied before for complex atoms or, at least, not within a relativistic framework. 
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TEST RUN OUTPUT 

** Example 1 ** 

Generation of radial Green's functions 

i E ( eV) j overall progress 

1 -l.OOOOOOOE+04 s 507. 

2 -1.5000000E+04 d- 100'/, 



Tests on the accuracy of the Green's functions by means of overlap and normalization integrals: 



i 


E 


( eV) 




Overlap integrals 
(Greens) 1 nj (Sal vat 95) > 


Normalization 
1 Inj (Greens) 1 1 


1 


-1, 


.OOOOOOOE+04 


Is 


1 . 000020E+00 


1 . 000057E+00 


1 


-1, 


,0000000E+04 


2s 


1.001925E+00 


1 . 003888E+00 


2 


-1, 


, 5000000E+04 


3d- 


9.997536E-01 


9.995404E-01 


2 


-1, 


.5000000E+04 


4d- 


1.005937E+00 


1 . 012284E+00 



Write the Green's functions to the .rgf file; 
2 radial Green's functions with 5065446 bytes ... 

XGREENS complete . . . 



** Example 2 ** 

Generation of radial Green's functions 



i 


E 


( eV) 


j 


overall 


1 


-1, 


,0000000E+03 


s 


16% 


2 


-1, 


, OOOOOOOE+04 


s 


33% 


3 


-1, 


, 5000000E+04 


P 


50% 


4 


-1, 


, OOOOOOOE+04 


P- 


66% 


5 


-1, 


.5000000E+04 


d- 


83% 


6 


-1. 


.OOOOOOOE+04 


d 


100% 



Tests on the accuracy of the Green's functions by means of overlap and normalization integrals: 



i 


E 


( eV) 


<nj 


Overlap integrals 
(Greens) I nj (Sal vat 95) > 


Normalization 
1 Inj (Greens) I I 


1 


-1 


.OOOOOOOE+03 


Is 


1 


.000267E+00 


1 


.000540E+00 


1 


-1 


.OOOOOOOE+03 


2s 


1 


.000318E+00 


1 


.000641E+00 


2 


-1 


.OOOOOOOE+04 


Is 


1 


.000003E+00 


1 


.OOOOlOE+00 


2 


-1 


.OOOOOOOE+04 


2s 


1 


.000138E+00 


1 


.000278E+00 


3 


-1 


.5000000E+04 


2p 


9 


.993563E-01 


9 


. 987343E-01 


3 


-1 


.5000000E+04 


3p 


9 


.999064E-01 


9, 


.998404E-01 


4 


-1 


.OOOOOOOE+04 


2p- 


1 


.000384E+00 


1, 


. 000770E+00 


4 


-1 


.OOOOOOOE+04 


3p- 


1 


.000216E+00 


1, 


. 000435E+00 


5 


-1 


.5000000E+04 


3d- 


1 


. 000046E+00 


1, 


. 000107E+00 


5 


-1 


.5000000E+04 


4d- 


1 


.001539E+00 


1 


.003181E+00 


6 


-1 


.OOOOOOOE+04 


3d 


1 


. 000084E+00 


1 


.000173E+00 


6 


-1 


.OOOOOOOE+04 


4d 


1 


.001030E+00 


1 


.002106E+00 



Write the Green's functions to the .rgf file; 
6 radial Green's functions with 59754187 bytes ... 

XGREENS complete . . . 
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Abstract 

Prom perturbation theory. Green's functions are known for providing a simple and 
convenient access to the (complete) spectrum of atoms and ions. Having these functions 
available, they may help carry out perturbation expansions to any order beyond the first 
one. For most realistic potentials, however, the Green's functions need to be calculated 
numerically since an analytic form is known onl\" for /'/■(;(; elect runs or for I heir mot. ion in 
a pure Coulomb held. Therefore, in order to facilitate Ihe use of (Troon's functions also 
for atoms and ions other than the hydrogen like ions, here we provide an extension to 
the Ratip program which supports the comput.ation of relativistic {one-electron) Green's 
functions in an — arbitrarily given — central-held potential V(r). Dilferent comput ational 
modes have been implemented to dehne these effective potentials and to generate t he radial 
Green's functions for all bound-state energies E < i). In addition, care has been taken 
to provide a user-friendly component of the Ratip package by utihzing features of the 
Fortran 90/95 standard such as data structures, allocatable arrays, or a module-oriented 
design. 
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PROGRAM SUMMARY 



Title of program: Xgreens. 
Catalogue number: To be assigned. 

Program obtainable from: CPC Program Library, Queen's University of Belfast, N. Ireland. 
Users may also down-load a tar-file of the program ratipOO.tar from our home page at the 
University of Kassel {http:/ /www. physik.uni-kassel.de/fritzsche/programs. html). 

Licensing provisions: None. 

Computer for which the new version has been tested: PC Pentium II, III, IV, Athlon. 
Installations: University of Kassel (Germany). 
Operating systems: SuSE Linux 8.2, SuSE Linux 9.0. 

Program language used in the new version: ANSI standard Fortran 90/95. 

Memory required to execute with typical data: On a standard grid (400 nodes), one central- 
field Green's function requires about 50 kBytes in RAM while approximately 3 MBytes are 
needed if saved as two-dimensional array on some external disc space. 

No. of bits in a word: Real variables of double- and quad-precision are used. 

Peripheral used: Disk for input/output. 

CPU time required to execute test data: 2 minutes on a 450 MHz Pentium III processor. 
Distribution format: compressed tar file. 

Keywords: Central— field Green's function, confluent hypergeometric function. Coulomb 
Green's function, Kummer function, multi-conhguration Dirac— Fock, regular and irregular 
solutions. 

Noiure of the physical probh-rn: In aloinic pc^'turbation theory, Green's functions may help 
carry out the summation over the complete spectrum of atom and ions, including the (sum- 
mation over the) bound states as well as an integration over the continuum [1]. Analytically, 
however, these functions are known only for free electrons (V(r) = 0) and for electrons in a 
pure Coulomb field (V(r) = —Z/r). For all other choices of the potential, in contrast, the 
Green's functions must be determined numerically. 

Method of solution: Relativistic Green's functions are generated for an arbitrary central-field 
potential V(r) = — Z(r)/r by using a piecewise linear approximation of the effective nuclear 

charge function Z(r) on some grid ri [i — 1 A^): Z;(/') = Zoi + Z\_,r. Then, following 

McGuire's algorithm [2], the radial Green's functions are constructed from the (two) linear- 
independent solutions of the homogeneous equation [3] . In the computation of these radial 
functions, the Kummer and Tricomi functions [4] are used extensively. 

Restrictions onto the complexity of the problem: The main restrictions of the program concern 
the shape of the effective nuclear charge Z(r) = — rV(r), i. e. the choice of the potential. 
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and the allowed energies. Apart from obeying the proper boundary conditions for a point-like 

nucleus, namely, Z(r — > 0) = .^rmc > and Z(r — * oo) = Z^^^c — A'eiectrons 0, the first 
derivative of the charge function Z(r) must be smaller than the (absolute value of the) energy 
of the Green's function, ^ < \E\. 

Unusual features of the program: XgreeNS has been designed as a part of the Ratip pacfcige 
[5] for the calculation of relativistic atomic transition and ionization properties. In a short 
dialog at the beginning of the execution, the user can specify the choice of the potential 
as weU as the energies and the symmetries of the radial Green's functions to be calculated. 
Apart from central-field Green's functions, of course, the Coulomb Green's function [6] can 
also be computed by selecting a constant nuclear charge Z(r) = Z^j^. In order to test the 
generated Green's functions, moreover, we compare the two lowest bound-state orbitals which 
are calculated from the Green's functions with those as generated separately for the given 
potential. Like the other components of the Ratip package, Xgreens makes careful use of 
the Fortran 90/95 standard. 

References: 

[1] R. A. Swainson and G. W. F. Drake, .1. Phys. A 24 (1991) 95. 
[2] E. J. McGuire, Phys. Rev. A 23 (1981) 186. 

[3] P. Morse and H. Feshbach, Methods of Theoretical Physics (McGraw-Hill, New York 1953), 
Part. 1, p. 825. 

[4] J. Spanier and B. Keith, An Atlas of Functions (Springer, New York, 1987). 

[5] S. Fritzsche, J. Elec. Spec. Rel. Phen. 114-116 (2001) 1155. 

[6] P. Koval and S. Fritzsche, Comput. Phys. Commun. 152 (2003) 191. 
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LONG WRITE-UP 



1 Introduction 

The use of Green's function has a long tradition for solving physical problems, both in classical 
and quantum physics. From perturbation theory, for instance, these functions are known for 
providing a rather simple access to the complete spectrum of a quantum system and, hence, to 
facilitate the computation of perturbation expansions beyond the first order in perturbation 
theory. Applications of the Green's functions can be therefore found not only in atomic and 
molecular physics but also in quantum optics, field theory, sohd-state physics, and at various 
place elsewhere. 

In atomic physics, however, the use of Green's function methods was so far mainly restricted 
to describe the motion of electrons in a pure Coulomb field, i. e. to the theory of hydrogen and 
hydrogen-like ions. For these one-electron systems, calculations have been carried out, for 
example, for the two-photon [1, 2] and multi-photon ionization [3, 4], the two-photon decay 
[5], the second-order contributions to the atomic polarizabihties [6] as well as for determining 
radiative corrections [7, 8]. 

Less attention, in contrast, has been paid to utilize Green's functions for non-Coulomb 
fields or for describing the properties of many-electron atoms and ions. Unlike to the gen- 
eration of the — bound and free-electron — wave functions to the Schrodinger and Dirac 
equation, for which a number of programs are available now also within the CPC-library 
[9, 10], there is almost no code freely available which helps generate the relativistic central- 
field Green's functions. As known from the literature, however, nonrelativistic central-field 
Green's functions were constructed by McGuire [11] and by HuiUier and coworkers [12], and 
were successfully applied for studying the multi-photon ionization of valence-sheU electrons 
in alkali atoms. Therefore, in order to facilitate the use of relativistic central-field Green's 
functions for atomic computations, here we describe and provide an extension to the Ratip 
package which calculates these functions (as the solution of the Dirac equation with a (5— like 
inhomogeneity) for an arbitrary central field V(r). 

In the following section, we start with summarizing the basic formulas for the computation 
of relativistic central-field Green's functions. Apart from a brief discussion of the Dirac Hamil- 
tonian, this includes the defining equation for central-field Green's function and the separation 
of the three-dimensional Green's function into radial and angular parts. However, since the 
separation has been discussed in detail elsewhere in the literature [13], we restrict ourselves to 
a short account on that topic and mainly focus on the computation of the radial components 
of the Green's functions. Section 3 later describes the program structure of Xgreens, its 
interactive control and how the code is distributed. Because the Xgreens program is de- 
signed as part of the RATI? package, we have used and modified several modules which were 
published before along with other components of the program. Section 4 explains and displays 
two examples of Xgreens, including (a) a dialog in order to calculate a central-field potential 
from Grasp92 wave functions [9] and (b) the generation of the radial Green's functions if the 
potential is loaded from an external file. Finally, a short summary and outlook is given in 
section 5. 
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2 Theoretical background 



2.1 Relativistic centred— field Green's functions 

Most naturally, the relativistic central-field Green's function can be considered as a general- 
ization of the (relativistic) Coulomb Green's function if, in the Dirac Hamiltonian^ 



-ica V + jit? + V{r) , 



(1) 



the Coulomb potential is replaced by some (arbitrarily given) central-field potential, V (r) = 
— f — > — As in the nonrelativistic case where the Green's function obeys the defining 
equation (H ~ E) Gir(r. r') = S(r — r') , the relativistic Green's function is given by a 4 x 4 
matrix [14] which satisfies the inhomogeneous equation 



{Ho(r)-E-c') GE(r,r') 



5(r-r')l4 



(2) 



with I4 being the 4x4 unit-matrix and where, as usual in atomic str\icture theory, E refers 
to the total energy of the electron but without its rest energy c^. Moreover, since in polar 
coordinates a central-field potential Z(r) in the Hamiltonian (1) does not afi^ect the separation 
of the variables, the central-field Green's function has the same radial-angular representation 
as in the pure Coulomb case [13] 



, „ ^ I ( Ekl:(r,r')n.^(v)ni^(i') -ig|f(r,r')n.™(?)nL.„(n\ _ 

Gi5(r,r') = > ^ . (3) 

^'''^ Vig|J('-,'-')J^-cm(f)nL(f') g|f('-,'-')J^-cm(f)nl.„(f') / 



where Q,f^jn[v) = ./^) denotes a standard spherical Dirac— spinor and k. — ± ( j + 1/2) 

for i = J ± 1/2 is the relativistic angular momentum quantum number; this number carries 
information about both the total angular momentum j as well as the parity (—1)' of the Green's 
function. While the summation over k = ±1, ±2, ... runs over all (non-zero) integers, the 
summation over the magnetic quantum number m = —j, + is restricted by the 

corresponding total angular momentum. Moreover, the radial part of the central-field Green's 



function 

Vgll ('■.'•') g|f ('•,?■') 

satisfies the equation 



in (3) can be treated simply as a 2 x 2 matrix function which 



/[ "Z(r) 



V 



aE 



K d 
r dr 
aZ(r) 



- aE 



a5{r — r') I2 , 



(4) 



with I2 now being the 2x2 unit— matrix. Note that in Eq. (4), a refers to Sommerfeld's fine 
structure constant and that, in order to keep the equations similar to those for the Coulomb 
Green's functions [13], we make use of a nuclear charge function Z(r) = — rV(r) to define the 
central-field potential, instead of V(r) explicitly. In the radial-angular representation (3), two 
superscripts T and T' were introduced to denote the individual components in the 2x2 radial 
Green's function matrix. These superscripts take the values T = L or T = 5 to refer to 



^Here and in the following, we use atomic units (me = h = e^/ineo = 1) if not stated otherwise. 
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either the large or small component, respectively, when multiplied with a corresponding (2- 
spinor) radial solution of the Dirac Hamiltonian (1). For a pure Coulomb potential, Z{r) = 
Zgg, an explicit representation of the (four) components ggj (r, r') of the radial Green's 
function can be found in Refs. [13, 14, 15]. 



2.2 Generation of radial Green's functions 



When compared with the Coulomb Green's functions, not much need to be changed for the 
central-field functions in Eqs. (3) and (4) except that the nuclear charge Z = Z(r) now depends 
on r and, hence, that analytic solutions to these equations are no longer available. Therefore, 
to find a numerical solution to Eq. (4), let us first mention that this matrix equation just 
describes coupled equations for the two independent pairs (g^^, gj^) and (g^^, gg^) of the 
radial components. For example, if we consider the first pair (gj;^, gj^) of radial components 
of the Green's function, it has to satisfy the two equations 



Q Z(r) 



■ aE 



dr 



&EAr, r') + 



aZ(r) 



_ 1^ 
dr 

-aE 



gl^(r,r') 
Slk{r,r') 



aS{r — r'), 



0, 



(5) 



(6) 



and a similar set of equations holds for the second parr (gjf , ggf )• In the following, therefore, 
we will discuss the algorithm for solving Eqs. (5) and (6) but need not display the analogue 
formulas for the pair (gjf , ggf )• 

Inserting gj^ (r, r') from Eq. (6) into (5), we arrive at the second-order inhomogeneous 
differential equation 



d 



Sk!:ir,r') 



a Z(r) 
+ —+aE 



a Z(r) 



+ gii-(r,r') = aS{r-r') 



(7) 



for the component g^^ir, r'). Solution of this equation can be constructed as product of two 
hnearly independent solutions [14] 



g|^(r,r') = Mi;^{min(r, r')) ■ W|,^ {max(r. r')) 



(8) 



which is regular at the origin, and W^^ (r) a solution regular at infinity. Below, we will obtain 

these functions following the numerical procedure as suggested by McGuire [11]. 

For this, let us start with approximating the nuclear charge function Z(r) on some grid 
n, i = 1, . . . , imax in terms of a set of straight lines 



Zi(r) = Zoi + Zur, 
and from the homogeneous part of Eq. (7) 



for ri < r < ri+i , 



.\[dr + r) g^-^'') 



> 9r 



2 aZ,(r) 



+ aE 



aZi(r) 



(9) 



(10) 
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as given for the i-th interval [n, ri+i] of the grid. In this Eq., we can drop the second 
argument r' in the component g^^, since it now appears only as a parameter, and replace it 
by the superscript i to denote the particular piece of the grid for which we want the solution. 
From the approximation (9), moreover, we see that the nuclear charge function Z.,;(r) within 
the j— th interval gives rise to a pure Coulomb potential Zoi/r and a constant which is simply 
added to the energy: E ^ Z^i + E. We therefore find that, within the given interval, the 
regular and irregular solutions at the origin in (8) can both be written as Unear combinations 
of the corresponding solutions for a Coulomb field 

mL(0 = /i,iMi!°-'°"'^r) + /i,2WiC°'^'°-b(r) (11) 
WL(0 = Si,iM^,°-'°-''(r) + s,,2W^°-'°-^r) (12) 

with constants {/»,!, /i,2, 9i,i, Si, 2} which still need to be determined. For the Coulomb po- 
tential, the functions M^^'"^^°'°^ (r) and Wg^'"''''™''(r) are known analytically and given by 
[14] 

^j,Couiombj^j = r'' e"'"- [U M{-ti + 1, 2si + 1, 2qir) 

+ (k - Zoi/qi) M(-tj, 2si + 1, 2qir)] , (13) 

^i. Coulomb = r*' e-*'- [(« + Zoi/qi) V{-ti + 1, 2si + 1, 2qir) 

+ \J{-ti,2si + l,2qir)] , (14) 

where M(a. h, r) and U(a, h, r) denote the Kummer and Tricomi functions [16, 17], respec- 
tively, and where the quantities s^, ti and Qi are given by 



Vl-((B + Zii)a2 + 1)2 



(15) 

V-(£; + Zii)((B + Zu)a2 + 2) . (16) 



The set of constants /i2i 9i:v-, — 1, '-max} f^^n be determined from the fact 

that the two functions Mj^(r) and Wj:^(r) in ansatz (8) as well as their derivatives need 
to be continuous in r, and that they behave regularly at the origin or at infinity, respectively. 
The constraint of being regular at the origin, for instance, requires the coefficients /i=i,i = 1 
and /j=i,2 = and can be used together with the continuity M^^(r) and Mg^^(r) , 

kiM^^in) + fi,2W^^{n) = /i+i,iMi+i'C(ri) + /i+i,2W^i:''°(ri), (17) 

/i,i^if('-i) + /i,2<!;°(n) = /i+i.iMi,';^i'°(ri) + /i+i.2W'El+''°(ri). (18) 

[where the superscript C here refers to the Coulomb functions in Eqs. (13) and (14)] in order 
to determine all the coefficients fij up to a normalization constant. A similar recurrence 
procedure also applies to the coefficients gij, but by starting from 'infinity', that is with 

fi'imax,! — and .giniax,2 = 1, and by going backwards in the index i towards the origin. 

To determine finally the normalization of the radial component gj^(y, r') , e. g. of the 
coefficients {fij, 9i,j}-, we may return to Eq. (7) and re-write it in the form 

(ho - e) g|^(r,r') = S{r-r'). (19) 
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Taking the integral over latter equation, we see that 

f^\ho-E) g^^(r, r') dr = f 5{t - r') dr = 1 
Jr'—e Jr'—s 

should be valid for £ — » + and for any r' and, hence, that — by using Eq. (7) and carrying 
out some algebraic manipulations — the derivative 



+ (20) 



'jumps' at r = r' . We can use the right hand side of Eq. (20) together with the derivative 
of Eq. (8) to determine the normalization constants c/ and Cg with which the coefficients fi^j 
and Si J need to be multiplied in order to obtain the normalized radial component gg^ (r, r') . 

Having constructed gj;^ (r. r') piecewise as sohition to Eq. (7), we may obtain the second 
component gj'^ (r, r') of this pair simply from Eq. (6), 



d^Ar,r')= \ , (21) 

- + ^-^ + aE 

a r 

where the derivatives of the Kummer and Tricomi functions in expressions (13) and (14) can 
be calculated by means of standard formulae [17] 

M'{a.,b,z) = ^M(a + l,6 + l,z), U'(a,6,z) = -aV{a, + l,b + 1, z) . 

In practise, both functions M{a,b,z) and lj{a,b,z) are required only for real arguments a, 6, 
and z but need to be calculated with a rather sophisticated algorithm [18] in order to ensure 
numerical stability and to provide sufficiently accurate results. A similar numerical procedure 
has to be carry out also for the second pair (gjf , ggf ) of radial components. This gives 
rise of course to another set of coefficients {fij, Qij} and, finally, to the full radial Green's 
function from Eq. (4). The 4 iniax coefficients for the (piecewise) regular and irregular solutions 
in ansatz (11-12) and (8) certainly provide — together with a few numerical procedures — 
the most compact representation of the radial central-held Green's functions. For the further 
computation of matrix elements and atomic properties, however, these radial components are 
usually represented (and stored) at some (2-dimensional) grid in r and r' as we will discuss 
below. 

2.3 Tests on the accuracy of the radial Green's functions 

To make further use of the Green's functions in applications, it is necessary to have a simple 

test on their numerical accuracy. For the radial-spherical representation of these functions 
as defined in Eq. (3), such a test is easily constructed since the Green's function contains the 
information about the complete spectrum of the Hamiltonian (1) and, hence, about any of its 
eigenstates V'n(r). Making use of the well-known expansion of the Green's function 

G.(r,r')^E^^#^ 
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(23) 



in terms of the eigenstates of the Hamiltonian, the relation 

^„(r') = [En -E) J dr 4(r) Geir, r') (22) 

can be easily derived from the orthogonality of the eigenfunctions and holds for all energies 
E ^ En- With the representation (3) in mind, the relation (22) can be written also in terms 
of the radial wave and Green's functions 

where P„k (r) and Q„k (r) are the large and small components of the (Dirac) radial 2-spinors. 
As indicated by the tilde on the left-hand-side of Eq. (23), therefore, this relation may serve 
as a test on the accuracy of the Green's functions if, for instance, the radial components from 
both side of the equation are compared with each other or if some proper overlap integral is 
calculated. 

To test on the precision of the Green's functions in Xgreens, we generate the radial 
bound-state wave function Pn/s(r) and Qn^ir) as solution of the radial Dirac equation 



■ aE 



dr 



with the same potential 



2 aZ(r) 



' dr 
-aE 



Q,Ur) 







(24) 
(25) 



_Z(r) 



as apphed before and by making use of the program by Salvat 
etal [10]. This solver has been embedded in our code and is utihzed in order to calculate the 

relativistic components. With these radial functions, we then compute the two integrals in 
(23) for obtaining P„^ (r) and Q„„ (r), respectively, as a function of r. As default in Xgreens, 
we use the overlap integral 



y ( Punir) P„^(r) + Q„„(r) Q„^r) ) dr 



(26) 



for the two lowest principal quantum numbers n of the (given) symmetry k, together with the 
corresponding normalization integral 



(P„„(r)P„„(r) + Q„„(r)Q„,.(r)) rfr, 



(27) 



as a numerical measure on the accuracy of the generated Green's functions. These integrals 
are displayed explicitly by the program (on demand). For a standard (logarithmic) grid with 
about 300 nodes, the overlap integral (26) is typically within the range 10~^ ... 10~^. For 
such a grid, the accuracy is limited in the computations by the linear interpolation of the wave 
and Green's functions. The accuracy can be increased however for a larger number «max of 
grid points as shown in section 4. 



3 Program structure 
3.1 The RATIP package 

Similar as Grasp92 [9] was designed for generating the wave functions within the multi- 
configuration Dirac-Fock (MCDF) model, the Ratip package [19, 20, 21] is organized as a 



9 



suite of program components to calculate a variety of (relativistic) atomic transition and 
ionization properties. Among other features, the components of RATI? support for instance 
investigations on the autoionization of atoms and ions, the interaction w4tli the radiation 
field, the parametrization of angular distributions in the emission of (electrons and photons, 
or the analysis of interference effects between radiative and non-radiative processes. In quite 
different case studies, Ratip has helped analyze and interprete a large number of spectra and 
experiments. In order to provide efficient tools for atomic computations, it also incorporates 
a number of further components (beside of the main program components for calculating 
certain physical properties) which facilitate the transformation of wave functions between 
different coupling schemes or the generation of continuum orbitals and angular coefficients. 
With the development of the XOKEENS program, we now provide an additional component to 
generate the relativistic central-field Green's functions within the Ratip environment. This 
is rather independent of their later 'use' where different properties such as the two-photon 
ionization and decay or various polarizabilities of atoms and ions might be calculated with 
the help of these functions. Owing to the large number of possible applications, however, a 
support of certain properties will largely depends on the requests by the users and on our 
further experience with the code. 

Not much more need to be said about Ratip's overall structure. A detailed account on its 
present capabilities has been given previously [22]. \Mtli regard to the further development 
of Ratip, we just note that our main concern now pertains to a long life-cycle of the code 
and to an object-oriented design within the framework of Fortran 90/95. With the present 
set-up of the Xgreens component, we continue our effort for providing an atomic code which 
is prepared for future applications in dealing with open-shell atoms and ions. 

3.2 Data structures and program execution 

Certainly, the major purpose of Xgreens is the computation of (one-electron) radial Green's 
functions for some — specified or externally given — central-field potential, which can later 
be utilized for calculating radial matrix elements of the type 

UZI*''" A; E; k,A) = J £ gl'ir) iK(kr) gif (r, r') y^{kr') g^"(r') dr dr' , (28) 

where jA(^^) denotes a spherical Bessel function, gj^^ (r, r') one of the components of the 
radial Green's function in (4), and where o and d refer to the radial wave functions of some 
bound or free-electron state, respectively. Matrix elements of this type play a key role, for 
instance, for studying two-photon ionization (bound-free) or decay (bound-bound) processes. 
Since the (four) components of the radial Green's functions depend on two radial coordinates, 
r and r', a rather large amount of data need usually to be calculated by means of Xgreens. 
These data are finally stored in an external (.rgf) radial Green's function file as will be 
discussed in the next subsection. 

To generate and utilize the Green's functions, two derived data types TGreens_single_rgf 
and TGreensjrgf have been introduced as shown in Figure 1. Using these data types, a single 
Green's function is kept internally in terms of its expansion coefficients from Eqs. (11) and 
(12) for the large-large and a similar set of coefficients for the small-small component on the 
given grid. Apart from the particular representation of the grid, these data structures contain 
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type : : TGrGGns_single_rgf 

! Data structure type to keep a single Green's function 

real (dp) :: energy 

integer : : kappa 

! Coefficients for the regular and irregular solutions 

real (dp), dimensionC : , : ) , allocatable :: fLLl, fLL2, fSSl, fSS2 

! Values of the regular and irregular solutions and their derivatives 
realCdp), dimensionC:), allocatable ;: mL, wL, mS, wS, & 

mLp, wLp, mSp, wSp 

end type TGreGns_single_rgf 
( 

! 

type : : TGreens_rgf 

! Data structure type to save the Green's functions 

r 

typeCTGrGens_singlG_rgf ) , dimensionC:), allocatable :: gf 

integer :: interpolation_mode , mtp 

realCdp), dimensionC:), allocatable :: yp, zO, zl 
end type TGreens.rgf 

Figure 1: The derived data structures TGreens-single jrgf and TGreensjrgf to generate, store, and 
interpolate a central-field Green's function. 

the energy and symmetry of the radial Green's function, the nuclear charge function(s) as well 
as information about the mode of interpolation of the radial components between the grid 
points and the maximum tabulation point mtp. In Xgreens, a variable of type(TGreens_rgf) 

is used in order to store the information about all the requested functions, and a few additional 
procedures are provided to calculate from these structures the values of the radial components 
for any set of arguments. 

The basic steps in the execution of the program are very similar as for other components 
of the Ratip package. At the beginning, all the necessary input data are read in by an 
(intoracti\'(;) dialog. The requested Green's functions are then calculated in turn and tested 
for their accuracy. The (main) output of the program is a formatted Ascil file which provides 
an interface for further apphcations, either within the Ratip environment or as worked out 
by the user. 

3.3 Interactive control and output of the program 

Like the other components of the Ratip package, Xgreens is controlled by a dialog at the 
beginning of the execution. In this dialog, all information need to be specified for obtaining the 
central-field potential as well as about the number and symmetry of the radial Green's func- 
tions to be generated. Owing to various choices for defining the spherical potential, however, 
slightly different dialogs may occur during the execution. An example is shown in Figure 2, 
where first a Hartree-plus-statistical-exchange (HX) potential due to Cowan [23, 24] is gen- 
erated from the ground-state wave functions of atomic gold, and then two radial Green's 
functions are calculated within this potential. 

Apart from the nuclear charge, given in terms of a Grasp92 .iso isotope data file, the 
dialog first prompts for specifying (or confirming) some basic parameters such as the speed of 
light, the grid parameters as well as for the energy units, in which the further input and output 
is done by the program. This is followed by a prompt for determining the central-field potential 
for the radial Green's functions. At present, Xgreens supports four models including a pure 
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Coulomb potential (for the calculation of Coulomb Green's functions) as well two potentials 
(HFS and HX) to incorporate also the exchange interaction in some approximate form. In the 
Hartree model, in contrast, only the direct potential of some atomic level is applied in the 
later computation of the Green's functions. In all these cases, however, the wave functions 
of the selected state have to be specified in terms of the corresponding Grasp92 files for 
the configm-ation state functions (CSF), mixing coefficients as well as the radial orbitals. If 
computed internally, the potential can be written also to disc and used in further applications. 
Beside of the internal set-up of the potential, in addition, a user-defined potential can be 
given also exphcitly and read in by the program. 

Although the nuclear charge is specified by means of a . iso file (in line with all other 
components of Ratip), the Green's functions are always generated for a point-like nucleus in 
order to 'utilize' their regular behaviour of the radial components at the origin [cf. section 2.2]. 
For each radial Green's function, then the energy and the one-particle symmetry k need to 
be specified. The functions are finally written to the .rgf radial Green's function file whose 
name has also to be given at the end of the dialog. In this file, the radial Green's functions 
are given at the grid (in r and r') as specified originally. Moreover, to test the accuracy of the 
Green's functions, a number of overlap and normahzation integrals are calculated on request 
in order to 'compare' a few low-lying orbitals, as obtained from the Green's functions by an 
integration over r' , with independcnitU- g(;nerated functions from Salvat's program [10]. In 
average, the computation of a single Green's function requires about 2 min on a 450 MHz 
Pentium PC and approximately the same time for the 'test' integrals. In the latter case, the 
rather large demand on CPU time arises mainly from the 2zniax radial integrals in r'— space, 
which need first to be performed in order to obtain the radial orbital functions from the given 
Green's function components. 

3.4 Distribution and installation 

A program of Ratip's size cannot be maintained over a longer period without that certain 
changes and the adaptation of the code to recent developments (in either the hardware or 
the operating systems) become necessary from time to time. For the distribution of the code, 
therefore, we follow our previous style in that the Ratip package is provided as whole. The 
main emphasis with the present extension of the program, however, is placed on the design 
and implementation of the two (new) modules rabs_greens.f90 and rabs_utilities_2.f 90 
which contain the code for the generation and the test of the radial Green's functions. In 
addition, the four modules (rabs_special_f unctions . f 90, rabs_error_control . f 90, 
rabs_greens_external . f 90, and rabs_cpc_salvat_1995 . f ) became necessary and had to be 
appended to the code in order to support an accurate computation of the hyper geometric and 
a few related functions from mathematical physics as well as for providing Salvat's solver [10] 
for the generation of the radial wave functions (P(r), Q(r)) within a given potential. 

In the present version, the Ratip package now contains the source code for the 8 components 
Anco, Cesd, Lsj, Greens, Rcfp, Relci, Reos as well as the Utilities for performing a 
number of small but frequently occurring tasks. Together with the six additional modules from 
above, the overall program therefore comprises about 40.000 lines of code, separated into 22 
modules (apart from the main program components and three libraries). AU these source files 
are provided in the Ratip root directory. As before, there is one makefile for each individual 
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component from which the corresponding executable can be obtained simply by typing the 
command make -f make-component, that is make -f make-greens in the present case. For 
most components of the Ratip package, in addition, we provide a test suite in a subdirectory 
test-component of the root wliere component refers to tlie names above. For example, the 
directory test -greens contains all the necessary files to run the sample calculations which we 
will discuss below. 

Before, however, the makefiles can be 'utilized', a number of global variables need to be 
specified for the compilation and linkage of the program. In a Linux or Unix environment, 
this is achieved by modifying (and sourcemg) the script file make-environment which saves 
the user from adopting each makefile independently. In fact, the script make-environment 
just contains a very few lines for specifying the local compiler, the options for the compiler, 
as well as the local paths for the libraries. Under Windows, in contrast, not much help need 
to be given to the user since most compiler nowadays provide the feature to define 'projects'. 
In this case, it is recommended to read off the required modules and source files from the 
makefile of the program component and to 'declare' them directly to the project as associated 
with the given component. 

In the past year, the XCREE.NS program has been applied under several LiNUX systems 
using the Lahey Fortran 95 compiler. The Ratip program as a whole has been found also 
portable rather easily to other platforms such as IBM RS/6000, SuN OS, or to the PC world. 
The file Read. me file in the ratip root directory contains further details for the installation. 
Overall, however, it is expected that it will not be difficult to compile the program also under 
other operating systems. 

4 Test calculations 

To illustrate the use of the Xgreens program, we briefly discuss and display two examples. 
They refer to the generation of the radial Green's function for a few selected (one-electron) 
symmetries in atomic gold {Z = 79) and demonstrate how the accuracy of these functions 
can be increased by a proper choice of the radial grid. However, before the radial components 
can be generated, we need first to specify the potential as shown in Figure 2. For the present 
examples, we have chosen a HX (Hartree-plus-statistical-exchange) due to the work of Cowan 
[23, 24], starting from the ground-state wave functions of neutral gold. As usual for Grasp92 
and Ratip, these wave functions have to be provided in terms of the . csl configuration state 
list and the .mix configuration mixing files as well as the . out radial orbital file. All of these 
input files are provided also in the test -greens/ subdirectory of the Ratip root directory 
and, thus, can be used for the present tests. In Figure 2 moreover, the potential is saved to 
the central-field potential file z79-au-hx . pot and later re-utilized in Figure 3. The output of 
the two examples are shown in the Test Run Output below. 

Like for the wave functions, the 'quality' of the generated Green's functions becomes only 
fully apparent, if they are used for calculating observaUes which can be compared to experi- 
ment. Although calculation of observables is not the aim of the present work, it is still possible 
to test the accuracy of the generated Green's functions by re-calculating the radial orbitals (as 
discussed in section 2.3) and by comparing these orbitals with those as obtained from the inte- 
gration of the Dirac equation for the given central field. This 'comparison' is done in Xgreens 
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XGREENS: Calculation of relativistic central-field Green's functions 
with energies E < (Fortran 95 version) ; 
CO Copyright by P Koval and S Fritzsche, Kassel (2004). 

Enter a file name for the greens. sum file: 
z79-au-pot . sum 

Enter the name of the isotope data file: 
z79. iso 

loading isotope data file . . . 

... load complete; 
The physical speed of light in atomic units is 137.0359895000000 ; 

revise this value ? 

The parameters of the radial grid 1 (first grid) are 
mt = 2.i77968408335618E-04 ; 
h = 6.250000000000000E-02 ; 
hp = O.OOOOOOOOOOOOOOOE+00 ; 
n = 390 ; 
revise these values ? 

Which units are to be used to enter and to print the energies: 
A : Angstrom; 

eV : electron volts; 

Hartree : Hartree atomic units; 
Hz : Hertz; 

Kayser : [cm**(-l)] ; 

eV 

Determine type of central-field potential to generate: 
Load : Load a potential from a .pot file 
Coulomb : (pure) Coulomb potential 
Hartree : Direct potential for a given ASF 
HFS : Hart ree-Fock- Slater method 

HX : Hartree-plua-statistical-eXchange (Covan 1967) 

HX 

Enter the name of the .csl (configuration states list) file: 
au.csl 

Loading configuration symmetry list file . . . 

There are 22 relativistic subshells; 

there are 3 relativistic CSFs; 

... load complete. 
Enter the name of the GRASP92 mixing coefficient file: 
z79-au . mix 
Loading mixing coefficients file . . . 

. . . load complete; 
Enter the name of the ruf data file: 
z79-au.rwf 

. . . load complete; 

Summary of all ASF to facilitate the selection of proper bound-state levels 
i Level J"P Energy (a.u.) Energy (eV) 



1 i 1/2 + -1.90397E+04 -5.18096E+05 
Enter ASF level number for which the potential should be calculated: 

1 

Enter the subshell (e.g. Is, 2p-, ...) which is specific to Couan's HX method: 

Is 

Write the generated potential to disc ? 

y 

Enter the name of .pot central-field potential file: 
z79-au-hx . pot 

Radial Green's functions will be calculated within the range <= r, r' <= 4.797 a.u. 
with the boundaries of the nuclear charge function ZCO) = 79.00 and Z(r_max) = 1.951 

In the given potential the one-electron energies (in eV) are: 
Is: -8.131E+04 

2s: -1.453E+04 2p-: -1.402E+04 2p: -1.214E+04 

3s: -3.479E+03 3p-: -3.232E+03 3p: -2.814E+03 3d-: -2.382E+03 3d: -2.289E+03 

Enter in turn the energies E (in eV) and symmetries of the radial Green's functions to be generated; 
Enter another energy emd symmetry (<cr> if done) : 
-10000 s 

Enter another energy and symmetry (<cr> if done) : 
-15000 d- 

Enter another energy and symmetry (<cr> if done) : 

Check the accuracy of the Green's functions after generation ? 

y 

Enter the name of .rgf file to save the Green's function: 
z79-au-pot .rgf 



Figure 2: Interactive dialog for calculating the effective potential and for generating the radial Green's 
functions. 



14 



The parameters of the radial grid 1 (first grid) are 
rnt = 2.1779684083356i8E-04 
h = 6.250000000000000E-02 
hp = O.OOOOOOOOOOOOOOOE+00 
n = 390 ; 
revise these values ? 

y 

Enter rnt : 
2 . 1779684083356 18E-04 

Enter h: 
3.i25E-2 

enter hp: 



enter n: 
780 

Which units are to be used to enter and to print the energies: 
A : Angstrom; 

eV : electron volts; 

Hartree : Hartree atomic units; 
Hz : Hertz ; 

Kayser : [ciii**C-l)]; 

eV 

Determine type of central-field potential to generate: 
Load : Load a potential from a -pot file 
Coulomb : (pure) Coulomb potential 
Hartree : Direct potential for a given ASF 
HPS : Hartree-Fock- Slater method 

HX : Hartree-plus-statistical-eXchange (Cowan 1967) 

Load 

Enter the nEime of a .pot file with a central-field potential 
z79-au-hx.fld 

Radial Green's functions will be calculated within the range <= r, r' <= 5.269 a.u. 
with the boundaries of the nuclear charge function Z(0) = 79.00 and Z(r_max) = 1.000 
In the given potential the one-electron energies (in eV) are: 
Is: -8.i31E+04 

2s : -1 .453E+04 2p- : -1 . 402E+04 2p: -1 . 214E+04 

3s: -3.479E+03 3p-: -3.232E+03 3p: -2.814E+03 3d-: -2.382E+03 3d: -2.290E+03 

Enter in turn the energies E (in eV) and symmetries of the radial Green's functions to be generated; 
Enter another energy and symmetry (<cr> if done) : 
-1000 s 

Enter another energy and symmetry (<cr> if done) : 
-10000 s 

Enter another energy and symmetry (<cr> if done) : 
-15000 p 

Enter another energy and symmetry (<cr> if done) : 
-10000 p- 

Enter another energy and symmetry (<cr> if done) : 

-15000 d- 

Enter another energy and symmetry (<cr> if done) : 

-iOOOO d 

Enter ajiother energy and symmetry C<cr> if done) : 

Check the accuracy of the Green's functions after generation ? 

y 

Enter the name of .rgf file to save the Green's function: 
z79-au-rgf .rgf 



Figure 3: Dialog for calculating the radial Green's functions with an increased number of radial points. 
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on request by evaluating for the two lowest principal quantum numbers n the (radial) overlap 
integrals of the re-calculated functions with the solutions by Salvat's program. In addition, 
the normalization integrals are also displayed as seen from the Test Run Output. Except 
of a few cases, the deviation of these integrals from 1.0 are typically well below 0.1 %. on the 
standard grid with about 300 mesh points in r and r'. For these deviations, there are two 
source of numerical 'inaccuracy' which arise from the test procedure and are not related to the 
generation of the radial Green's functions. They are caused by the fact that by using Salvat's 
solver [10], a third-order spline is used to represent the potential at intermediate points in r in 
contrast to the linear representation (9) for the generation of the Green's functions. Moreover, 
the integration of the radial integrals uses the standard Grasp92 procedure and, hence, is not 
adapted as well to Salvat's solutions. 

The accuracy of the Green's functions (and the test integrals) can be improved however 
if the number of grid is enlarged, i. e. the spacing between the mesh points reduced. In our 
second example, cf. Figure 3, we therefore adopt a grid with approximately twice the number 
of mesh points. Note that this requires four times as much storage as the radial components 
are functions in r and r'. In this test case, moreover, the radial Green's functions are generated 
for four different symmetries s, pi/2i Pzi2^ (^3/2) <^5/2- When the Test Run Output from 
this example is compared with those from example 1, we see that the accuracy in the overlap 
integrals in increased by about a factor of 5. For this example, the generation and test took 
about 30 minutes in total on a 450 MHz Pentium III. 

After the termination of Xgreens, the radial Green's functions file (e. g. z79-au-rgf . rgf 
from example 2) contains all the Green's function components which were generated during 
the execution. Since this is a formatted file, it can easily be read and manipulated by any 

text editor. This (radial Green's function) file stores the individual components gji'j (r, r') of 
the Green's functions as 2— dimensional arrays within a simple file structure. An (internal) file 
signature # DCFGF in the first line is followed by the mode of interpolation and the number of 
Green's functions which are provided by the file. Then, each Green's function is specified (in 
turn), starting with the energy and symmetry of the function and followed by a table of six 
rows where the first two rows refer to the coordinates r and r' (in atomic units) and and the 
other rows to the four component functions g_E^('^, "r'), SEfC^j^')) sJk(^'^0 gjf(^)^') 
this particular order. Although this format is convenient for the later use of these functions, 
it may lead to rather large data files, especially if some larger number of Greens' functions 
need to be generated. In practise, however, not much problems are expected with this file 
structure because disc storage became cheap recently and because, for most applications, we 
expect the Green's functions to be generated on demand by making use of the corresponding 
modules instead of keeping them in external files. Typically, the radial wave and Green's 
functions occur as part of matrix elements and, thus, first require an additional integration 
(over r and/or r') before any observable quantity is obtained. 

The format of the radial Green's function (.rgf) file can be used easily in order to plot 
the various radial components. In Figure 4, therefore, we display the nuclear charge function 
Z(r) and the t wo radial Green's function components g^f (r|),r) and g£^{v[).r) as function 
of (second) argument r and taken for ro = 0.204076 a. u. The radial Green's functions are 
calculated in the ground-state potential of atomic gold, as apphed in example 2, for the energy 
E = 1000 eV. As seen from Figure 4, the large-large component gj^ is continuous at r = r' 
while the small-large component gJJ jumps owing to Eqs. (20) and (21) at this position 
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Z = Z(r) 
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Figure 4: Nuclear charge Z — Z(r) and the radial Green's function components g^^(ro,r) and 
gj^{r(j,r) as function of r and for j'n — 0.204076 a.u. The radial Green's functions with ,s— symmetry 
{k — —1) arc calculated in the ground-state potential of atomic gold for the energy E — 1000 eV. 



in space. Both components, moreover, behave 'regularly' at the origin as was constructed 
explicity by Eq. {8). 

5 Outlook. Applications of central— field Green's functions 

Apart from fundamental interest in having the central-field Green's functions available for 
'relativistic' electrons, these functions are also useful for the perturbative treatment of atoms 
and ions. As discussed in subsection 2.3, namely, the Green's function provide a simple access 
to the complete spectrum of a quantum system and, hence, can be utilized for carrying out 
the summation over all the (one-particle) states of the spectrum as required in second- and 
higher-order perturbation theory. However, since we only provide the one-electron Green's 
functions, they can be applied for just those processes where the 'perturbations' are described 
in terms of one-particle operators. Perhaps, the most-studied perturbation of this type is the 
interaction of atoms and ions with the radiation field, which — within a sufficiently strong 
field — may lead not only to atomic photoionization but (in second order) also to the two- 
photon ionization and decay of atoms and ions and to several other processes. During the last 
two years, we utilized the central-field Green's functions mainly for exploring the two-photon 
ionization for the helium-like ions and for the inner-shell electrons from atomic neon and 
argon. 

Apart from the one-particle Green's functions (as discussed in this work), there are a 
number of further processes such as the double Auger decay, for which the two-particle Green's 
functions G£;(ri,r2; r'^^j^s) required in order to calculate the cross sections, decay rates 
and/or angular distributions. These two-particle functions are solutions to a defining equation 
similar to Eq. (2) but where the (one-particle) Dirac Hamiltonian Hd is replaced by the two— 



17 



particle operator Hy),i + -f^D,2 + V12, including the electron-electron interaction V12 explicitly. 
Even by making use of a proper radial-angular decomposition of such Green's functions, their 
radial part g^^ k2(''1' ''^2) would depend then on four radial variables and an overall rank 
fc, similar as known from the tensorial decomposition of the electron-electron interaction [25]. 
Therefore, an internal representation and generation of these (radial) functions appear rather 
infeasible in practise. An alternative to these two— particle Green's functions is given by the 
(so-called) modified pair functions ^ e .cki^l , ''^2) which satisfy the equation 

{Hd.i + Hu,2 + V12- *E.„{ri,r2) = W |*„(PJM)) , 

where Vt^ is a one- or two-particle operator which depends on the physical task to be solved 
and |\I/q(PJ^/)) is one of the (two-electron) solutions of the corresponding homogeneous 
equation. These modified pair functions have the advantage that they only depend on two 
radial variables similar as the one-particle Green's functions. On the other hand, of course, 
these functions now depend on some (initial) two-particle state |^'n(PJAf)) as well as on the 
particular choice of the interaction operator W and, hence, are less general when compared 
with the Green's function from above. We currently investigate the possibilities for generating 
also such modified pair functions within the framework of Ratip and for their (later) use in 
the description of the double Auger decay. But already the one— particle Green's functions 
from this work enables one to explore a number of atomic properties which have not been 
studied before for complex atoms or, at least, not within a relativistic framework. 
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TEST RUN OUTPUT 



** Example 1 ** 
*************** 

Generation of radial Green's functions 

i EC eV) j overall progress 

1 -1 . OOOOOOOE+04 s 507. 

2 -1 . 5000000E+04 d- 100°/. 

Testa on the accuracy of the Green's functions by means of overlap and normalization integral 



i 


E 


C eV) 


nj <nj 


Overlap integrals 
(Greens) 1 nj CSalvat 95)> 


Normalization 
1 Inj (Greens) 1 1 


1 


-1. 


, OOOOOOOE+04 


Is 


1.000020E+00 


1.000057E+00 


1 


-1. 


, OOOOOOOE+04 


2s 


1.001925E+00 


1.003888E+00 


2 


-1. 


.5000000E+04 


3d- 


9.997536E-01 


9.996404E-01 


2 


-1. 


.5000000E+04 


4d- 


1.005937E+00 


1.012284E+00 



Write the Green's functions to the .rgf file; 
2 radial Green's functions with 5065446 bytes ... 

XGREENS complete . . . 



*************** 
** Example 2 ■ 



of radial Green's functions 



i 


E ( eV) 


3 


overall 


1 


-l.OOOOOOOE+03 




ley. 


2 


-l.OOOOOOOE+04 




337. 


3 


-1.5000000E+04 


P 


soy. 


4 


-l.OOOOOOOE+04 


P- 


eey. 


6 


-1.5000000E+04 


d- 


832 


6 


-l.OOOOOOOE+04 


d 


looy. 



Tests on the accuracy of the Green's functions by means of overlap and normalization integral 



i 


E 


C eV) 


nj <nj 


Overlap integrals 
(Greens) ! nj (Salvat 95)> 


Normalization 
llnj (Greens) II 


1 


-1, 


,0000000E+03 


Is 


1 . 000267E+00 


1 . 000540E+00 


1 


-1, 


,0000000E+03 


2s 


1.000318E+00 


1.000641E+00 


2 


-1, 


, OOOOOOOE+04 


Is 


1 . 000003E+00 


l.OOOOlOE+00 


2 


-1, 


, OOOOOOOE+04 


2s 


1.000138E+00 


1 . 000278E+00 


3 


-1, 


,5000000E+04 


2p 


9.993563E-01 


9.987343E-01 


3 


-1, 


,5000000E+04 


3p 


9.999064E-01 


9.998404E-01 


4 


-1, 


, OOOOOOOE+04 


2p- 


1 . 000384E+00 


1 . 000770E+00 


4 


-1, 


.OOOOOOOE+04 


3p- 


1.000216E+00 


1 . 000435E+00 


5 


-1, 


.5000000E+04 


3d- 


1 . 000046E+00 


1.000107E+00 


5 


-1, 


. 5000000E+04 


4d- 


1.001539E+00 


1.003181E+00 


6 


-1. 


.OOOOOOOE+04 


3d 


1.000084E+00 


i.000173E+00 


6 


-1. 


, OOOOOOOE+04 


4d 


i.001030E+00 


i.002106E+00 



Write the Green's functions to the .rgf file; 
6 radial Green's functions with 59754187 bytes ... 

XGREENS complete . . . 
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